Lattice Dynamics and Melting of a Nonequilibrium Pattern 
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We present a new description of nonequilibrium square patterns as a harmonically coupled crystal 
lattice. In a vertically oscillating granular layer, different transverse normal modes of the granular 
square- lattice pattern are observed for different driving frequencies (fd) and accelerations. The 
amplitude of a mode can be further excited by either frequency modulation of fd or reduction of 
friction between the grains and the plate. When the mode amplitude becomes large, the lattice 
melts (disorders), in accord with the Lindemann criterion for melting in two-dimensions. 



Systems driven away from thermodynamic equilibrium 
often form patterns when forced beyond a critical thresh- 
old. Close to this bifurcation, the dynamics of the 
nonequilibrium patterns are well described by partial 
differential equations called amplitude equations, whose 
forms are universal [l|. However, these equations lose 
predictive power for larger forcing. As an alternative de- 
scription of patterns, Umbanhowar et al. |2| conjectured 
that interacting localized structures, "atoms" of the pat- 
terns, could be the building blocks of spatially extended 
patterns. In this Letter, we demonstrate that this specu- 
lation was correct: we describe and predict the behavior 
of a pattern in a particular system by treating it as a 
finite number of interacting elements. Such an approach 
can replace the description of patterns by partial differ- 
ential equations with a possibly simpler description based 
on a finite set of coupled ordinary differential equations. 

We use a well-known pattern forming system, a thin 
layer of oscillated granular material. We show that the 
dynamics (including all oscillatory motions and disor- 
dered states) of square patterns are analogous to the 
dynamics of a discrete lattice of harmonically coupled 
elements. The descriptive and predictive power of such a 
framework is especially useful for the patterns formed in 
granular materials because the governing equations are 
only now being rigorously tested p[ , and the exact forms 
of the amplitude equations are not yet established. 

Experiment. A layer of 0.17 mm diameter bronze 
spheres was oscillated vertically in an evacuated cell at 
driving frequency fd with non-dimensional peak plate ac- 
celeration T = A(2irfd) 2 /g, where A is the amplitude of 
the plate oscillation and g is the gravitational accelera- 
tion. For the layer depth studied (4 particle diameters), 
square patterns formed for 2.5 < Y < 4.0 and fd < 36 
Hz. The granular surface was imaged using low angle 
illumination that created bright regions at the peaks |4( . 
The scattered light was collected by a 256 x 256 pixel 
CCD camera. 

Lattice pattern. The square patterns, illustrated in 
Fig. da), oscillate subharmonically at /d/2; after each 
plate oscillation, a peak becomes a crater. At the phase 
in the plate oscillation cycle when the pattern ampli- 
tude is maximum, the pattern is composed of an array of 
peaks arranged in a square lattice connected by a network 
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FIG. 1: (a) A lattice pattern strobed at fa/2 and time av- 
eraged (r = 2.90 and fd = 25 Hz, parameters for which the 
lattice is almost motionless) . The image shows the entire con- 
tainer, 18 x 18 cm 2 , (b) Close up snapshot image (3x3 cm 2 ) 
at r = 2.90 and fd = 30 Hz, for which the lattice oscillates 
vigorously, (c) The evolution in time (with units T = 1/ fd) 
of the peaks in boxes A and B in (b) ; these peaks oscillate out 
of phase with a frequency about 10 times smaller than /d/2. 



of thin lines of particles. At maximum amplitude each 
peak typically contains several hundred particles. Im- 
ages are collected at this phase in the cycle. In the dark 
regions between the peaks, there are almost no grains. 
Thus, when strobed at /d/2, the pattern resembles a two- 
dimensional (2D) square crystal lattice made of discrete 
elements separated by lattice constant a. In this paper, 
we only consider the motion of the lattice pattern strobed 
at fd/2. 

Lattice oscillation. We will now show that the lattice 
analogy is more than superficial: the square patterns also 
exhibit the dynamics of a coupled lattice. In general, the 
center of mass of each peak (lattice element) oscillates 
around its equilibrium (time-averaged) lattice site. We 
find that for a range of control parameters, the oscillation 
is periodic; a particular example is illustrated in Fig.^Jb) 
and (c). The peaks in a row at an angle 7r/4 to the natural 
lattice direction maintain a constant separation of \pla 
as they oscillate. Such an oscillation of the pattern is 
analogous to a single excited transverse normal mode in 
the (1,1) direction of a 2D crystal lattice. We find that 
all periodic motions of the lattice resemble (1, 1)t modes 
(or the degenerate (1, — 1)t), with different wavevectors 
and oscillation frequencies. 

Dispersion relation. The evidence that the granular 
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FIG. 2: (a) Comparison of the measured (1, 1)t dispersion 
relation (o) at fd = 25 Hz and T = 2.75 with a lattice model 
(solid line). The wavenumber k is in units of 47t/(v / 2g^V), 
where a is the lattice constant and N = 34 is the number 
of (1,1) rows in our square container (N = 32 in Fig. QJ. 
Only N/2 positive k modes are shown since the lattice os- 
cillations are decomposed into periodic Fourier components. 
The dashed line denotes the edge of the first Brillouin zone 
(k = k gz) for the (1, 1)t modes, (b) and (c) show the mod- 
ulus of the spatial Fourier transforms, \I(k x , k y , /z,)|, at two 
lattice oscillation frequencies, fh — 2.3 Hz and fh — 1.2 Hz, 
respectively. The sidebands represent the spatial modulation 
of the lattice tt/4 from the basic square lattice direction — the 
(1, 1)t modes. The four peaks (•) that form the basic square 
lattice (found at fh = Hz) have a power about 50 times 
larger than the background. 

lattice behaves like a discrete coupled lattice is provided 
by the dispersion relation shown in Fig. Ufa). Since the 
lattice dynamics is in general more complicated than the 
periodic oscillation of a single mode shown in Fig. d 
the dynamical behavior is determined from the three- 
dimensional discrete Fourier transform of a time series of 
images, giving I(k x ,k y , fx/). Only certain values of os- 
cillation frequencies fx, and wavevectors (k x ,k y ) contain 
power (|/| 2 ), and these combinations yield the dispersion 
relation for the lattice in the (1, 1)t (or the degenerate 
(1, — 1)t) direction |5|. 

As a first approximation to the unknown 2D poten- 
tial between the lattice elements, following Kittel [(j, we 
assume that the modes can be decoupled in different crys- 
tal directions. Since we observe only the (1, ±1)t modes, 
we compare the data to a ID lattice model in which we 
assume that the transverse motion of each (1,1) row is 
harmonically coupled to its nearest (1,1) neighbors. As 
shown in Fig. Efa), the measured dispersion relation is 
well described by the dispersion relation of the ID lat- 
tice model with Jl = Jbz\ sin(/ca/(2v / 2)|, where fsz is 
the frequency at the edge of the Brillouin zone (propor- 
tional to the square root of the ratio of effective spring 
constant to the effective mass of a (1,1) row), a the lat- 
tice constant, and k is the magnitude of a wavevector in 
the (1,1) direction 7]. We have found that for a wide 
range of T and fd, the measured dispersion relations are 



fit well by the harmonically coupled lattice model; as T 
and fd are changed, fsz varies between 1.5 Hz and 2.5 
Hz, an order of magnitude smaller than fd. 

Excitation of normal modes. We have demonstrated 
that lattice dynamics describes the square patterns; 
therefore, we can view the periodic oscillations of the pat- 
tern as individually excited normal modes. In two ranges 
of T and fd, the oscillation of the pattern becomes very 
intense (Fig.JHJ). The two regions are characterized by ex- 
citation of specific normal modes: In region I, the power 
is dominant in the mode at the edge of the (1,1) Brillouin 
zone, k = ksz = 2^/ (V^a), while in region II, the power 
is dominant in a mode near the middle of the Brillouin 
zone. We do not know why these particular modes are 
excited. Away from the two regions of excitation, the 
lattice is nearly stationary, with small amplitude oscilla- 
tions of the lattice elements around the equilibrium sites. 
In these quiescent regions, the total power is small and 
is roughly independent of the mode number: the lattice 
acts as if it is in contact with a thermal bath. 
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FIG. 3: Normal modes of the lattice exhibit two regions of 
intense excitation, I and II. Top panel: the grayscale indi- 
cates the power in the dominant mode, which was obtained 
by integrating the spectral power above a noise background. 
Bottom panel: The wavevector of the dominant mode. The 
points are taken along the path in the top panel. 

Defect formation. In a real crystal lattice, if the am- 
plitude of the oscillations is large enough, defects form 
and melting occurs. We see the analogous effect in the 
square pattern lattice: near region II, the single mode 
oscillations are large enough to locally break the lattice 
and form defects, which either remain in the lattice or 
travel through the crystal and annihilate at the bound- 
ary. Thus, in region II the patterns do not display perfect 
long range order. 

The defect creation rate can be enhanced by resonantly 
exciting a specific normal mode, as shown in Fig. 0Ja). 
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FIG. 4: Defect creation and melting for three cases after a 
sudden change in the system parameters. The numbers are 
the number of oscillations after the change in conditions, and 
the insets show the structure factor (power spectrum) corre- 
sponding to the pattern, (a) At t = 0, frequency modulation 
with fmr = 2 Hz and / ms = 5 Hz was applied (T — 2.9, 
fd = 32 Hz), (b) At t = 0, the same frequency modulation as 
in (a) was applied for particles which had been cleaned and to 
which graphite had been added, (c) Molecular dynamics sim- 
ulation: at t = 0, the friction coefficient \i between the grains 
and the plate was reduced from 0.5 to (T = 3.0 and fd = 32 
Hz, for which the pattern weakly oscillates when \i — 0.5). 



We achieve this by a slow frequency modulation of the 
layer oscillation: the signal applied to the container has 
the form y = Asm(2,7rfdt + j 2 ^ sin 27r/ mr t), where fd is 
modulated at a rate f mr with a modulation span f ms . 
The lattice responds to the frequency modulation by os- 
cillating in a (1, 1)t normal mode with a frequency of 
exactly / mr /2, and the strength of the response increases 
as f ms increases. Several hundred oscillations after the 
modulation is turned on, the amplitude of the excited 
mode becomes quite large, and a few defects are created 
as the mode locally shears the lattice apart. 

Friction and melting. To further increase the rate of 
defect formation, we have reduced the friction between 
the particles and between the particles and the plate by 
adding approximately 0.2 grams of graphite powder to 
the 60 grams of freshly washed particles. This results in 
enhanced defect formation within the lattice near region 
II even without frequency modulation of the container || ; 
frequency modulation of these particles fully melts the 
pattern (see Fig. 0[b)), creating a disordered liquid-like 
pattern in which lattice elements are no longer confined 
to oscillate near their equilibrium positions. We hypoth- 
esize that the friction controls the mode amplitude by 
damping the in-plane sliding motion of the peaks against 
the container bottom. 

In the experiments we cannot vary the friction in a 
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FIG. 5: Melting occurs when the Lindemann ratio jm (see 
text) reaches approximately 0.1. Main figure: the evolution 
of 7m in simulations with V = 3.0 and fd — 32 Hz after the 
friction coefficient \i was changed at t = from \i — 0.5. 
The inset shows that the correlation length of the pattern for 
jji — reached the plateau region when jm ~ 0.1. 



controlled way, but this can be done in molecular dy- 
namics (MD) simulations. We previously found that an 
MD code for inelastic hard spheres yielded square pat- 
terns like those in Fig.Qfa) 9]. With a friction coefficient 
/jl = 0.5, the spatial patterns and their wavelengths ob- 
tained from MD simulations were in good accord with 
the experiments for a wide range of T and fd- We find 
that melting in the MD simulation occurs rapidly when 
the friction coefficient between the grains and the plate 
is reduced to zero, as Fig.^fc) illustrates. In this melting 
process the correlation length of the pattern £ decreases 
from about 5a to 2a (see inset of Fig. EJ). (£ was ob- 
tained by fitting an exponential to the envelope of the 
azimuthal average of the 2D autocorrelation function.) 
After reaching a minimum in about 30 plate oscillations, 
£ oscillates about 2a as defects continue to appear and 
disappear in the liquid-like pattern. Since the melting of 
the granular lattice is driven by the growth of the single 
excited mode, we cannot compare our observations with 
the well-known 2D melting theory of Halperin and Nel- 
son which requires uniform (thermal) heating of all 
modes. 

Lindemann melting criterion. In simulations of 2D 
crystalline lattices, melting has been found to occur when 
the Lindemann ratio, 7m = (\u m — u n \ 2 )/a 2 , exceeds 
0.1 [DIE^; here u denotes displacements of atoms from 
lattice sites, and the average is taken over all nearest 
neighbors m and n. The time evolution of 7m computed 
for the granular lattice for different values of \i is shown 
in Fig. We find that the lattice melts (£ reaches its 
plateau region) when 7m increases through the value of 
0.1, just as in 2D crystals. For \i = 0.5, the amplitude of 
the excited mode is too small for the lattice to melt, while 
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for fi = 0.1, the pattern oscillates with increasing ampli- 
tude until local melting events disrupt the long-range or- 
der (Fig. El — such behavior is analogous to the dynamics 
seen in the experiment near the peak of region II. Thus, 
the Lindemann ratio can be used as a predictive criterion 
for the loss of order in the square patterns. 

Conclusions. We have shown a close correspondence 
between the square patterns in a vibrated granular layer 
and a discrete set of harmonically coupled elements in a 
2D lattice. The applicability of the Lindemann criterion 
in the nonequilibrium system further confirms the util- 
ity of the lattice dynamics framework. In addition, we 
conjecture that the stability exhibited by the nonequilib- 
rium granular lattice pattern is an example of the "gen- 
eralized rigidity" typically found in equilibrium lattice 
systems and discussed by Anderson Whether such 
concepts and predictive criteria can be applied to other 
nonequilibrium systems is an open question. For exam- 
ple, it would be interesting to see if other nonequilibrium 
patterns that display excited modes [l4j exhibit a dis- 
persion relation and a melting scenario like that we have 
observed. Finally, we note that our observations of the 
effect of friction on melting of the lattice could help guide 
the development of kinetic and hydrodynamic theory of 
granular media, where the role of friction is not well un- 
derstood. 
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